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ABSTRACT 

In this letter we present the first 3D spectroscopic study of Ha emission in a z ~ 1 
field galaxy with an integral field unit. Using the CIRPASS spectrograph on Gemini- 
South we map the spatial and velocity distribution of Ha emission in the z = 0.819 
galaxy CFRS 22.1313. We detect two Ha emitting regions with a velocity separation of 
220 ± 10 km s^^. Combining the 2D map of Ha emission with HST F814W imaging, we 
determine a lower limit of 180 ± 20 km s~^ for the rotation velocity of this Afs(rest) ^ - 
21 galaxy. We note that our value is significantly higher than the rotation velocity 
of 120 ± 10km s~^ reported by Harden et al. (2003) for their long-slit spectroscopic 
study of this galaxy. Our lower limit on the rotation velocity is entirely consistent 
with no evolution of the rest B-band TuUy-Fisher relation. The position of this galaxy 
relative to the mean rest -B-band TuUy-Fisher relation of TuUy & Pierce (2000) is 
consistent with brightening of no more than ~lmag at z = 0.8. A larger integral 
field unit sample, without the uncertainties inherent to long-slit samples, is needed to 
accurately determine the evolution of the TuUy-Fisher relation out to z ~ 1. 

Key words: galaxies: evolution - galaxies: individual: CFRS22.1313 - instrumenta- 
tion: spectrographs - galaxies: spiral - galaxies: kinematics and dynamics 



1 INTRODUCTION 

The TuUy-Fisher (TP) relation describes the strong corre- 
lation between luminosity and maximum rotation velocity 
for disc galaxies (TuUy & Fisher 1977). This reflects a fun- 
damental relationship between the total mass of the galaxy 
and the mass contained in stars. The redshift evolution of 
this scaling relation provides a powerful test of galaxy forma- 
tion models. Hierarchical models of galaxy formation predict 
that the ratio of stellar mass to total mass would be simi- 
lar at all redshifts i.e. massive galaxies assembled through 
mergers of smaller galaxies. The 'classical' galaxy forma- 
tion models (e.g. monolithic collapse; Eggen, Lynden-Bell 
& Sandage 1962) predict a higher mass to light ratio at 
higher redshift, since the gas is still being converted to stars 
within the fully-formed dark matter halo. However, stellar 
population evolution might produce an offset in the oppo- 
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site sense (when based on observations at short rest-frame 
wavelength), since at high redshift a stellar population of 
identical mass is likely to be younger and therefore bluer 
than an equivalent population at low redshift. 

Investigation of the evolution of the rest frame _B-band 
TF relation at z > 0.5 has to date produced discrepant 
results with observations of only moderate luminosity evo- 
lution of <0.4 mag out to 2 ~ 1 (Vogt et al. 1996, 1997) 
compared with claims of much stronger evolution of ~ 2 
mag at z ~ 0.3 (Simard & Pritchet 1997; Rix et al. 1997). 
Previous high-z work has so far been limited to long-slit 
spectroscopy, which has inherent problems - missing light 
through slit losses and possible misalignment of the slit, po- 
tentially leading to erroneous results for the rotation velocity 
if not addressed properly. However, the relatively new tech- 
nique of 3D spectroscopy, used by Andersen & Bershady 
(2003) at low redshift, overcomes these difficulties. Using an 
integral field unit (IFU) we can efficiently obtain a com- 
plete census of the spatially extended line emission. We can 
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potentially map disc kinematics at z ~ 1 out to sufficiently 
large galactocentric radii to measure Vmax (without need for 
correction), as well as obtaining more detailed information 
on the spatial distribution of emission line regions. 

In this letter we present the first demonstration of near- 
IR integral field spectroscopy of a z ~ 1 field galaxy, using 
our new CIRPASS instrument (Cambridge IR PAnoramic 
Survey Spectrograph, Parry et al. 2000) on Gemini-South. 
By moving into the near-IR we were able to study the star 
formation in the z = 0.8 galaxy CFRS 22.1313 (Lilly et al. 
1995b) using the rest-optical Ha A 6563 A emission line, the 
same reliable tracer of star formation as used at low redshift. 
We investigate the disc kinematics through the Ha line and 
compare the results found here using an IFU to previous 
work using traditional long-slit spectroscopy. 

The layout of this letter is as follows. In section 2 we 
describe our target galaxy, the CIRPASS IFU observations 
and archival HST imaging. In section 3 we discuss the dis- 
tribution of star formation and the large-scale kinematics of 
the galaxy. Our conclusions are presented in section 4. We 
assume Ho = 70 km s~^Mpc~^, Qm = 0.3 and $1a = 0.7 
throughout this letter, unless otherwise stated. An angular 
scale of 1 arcsec corresponds to a physical distance of 7.57 
kpc ai z — 0.819 using the above cosmology. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 CIRPASS Spectroscopy 

Our target is the disc galaxy CFRS 22.1313 (part of the 
Canada-France Redshift Survey (CFRS) - Lilly et al. 1995a) 
and was specifically chosen for IFU study as it is spatially 
extended and has an accurately known redshift {z — 0.819), 
for which Ha appears between sky lines. It also has strong 
[OH] A 3727 A fiux, which should imply detectable Ha emis- 
sion (Kennicutt 1992), and benefits from deep HST imaging 
from which the disc structural parameters may be derived. 
Near- infrared integral field spectroscopy of CFRS 22.1313 
was obtained with CIRPASS on Gemini-South on the night 
of 2002 August 13. CIRPASS is a fibre-fed spectrograph with 
a 490 lenslet array with a variable lenslet scale and operates 
in the J- and iif-band (1-1.67 /im). Using the 0.25 arcsec di- 
ameter lenslet scale for these observations the array covered 
an area of 9.3 by 2.9 arcsec. The detector is a Ik x Ik Hawaii-I 
HgCdTe Rockwell array. The CIRPASS observations com- 
prised eight separate 1800 s integrations with a seeing of 
about 0.4 arcsec. An eight-point dither pattern was used to 
facilitate sky subtraction and the removal of bad pixels. The 
spectra were taken using a 400 1 mm^^ grating with resolv- 
ing power R — ~ 3000 and dispersion 2.28 A pix~^. The 
wavelength coverage was 1-1.24 /im, targeting the redshifted 
Ha emission line, which appeared in a region free of OH sky 
emission lines. 

The data were reduced using the cirpass iraf package^ . 
For each 1800 s integration the detector was read out non- 
destructively 20 times at the start of the integration and 
at 600 s, 1200 s and 1800 s. The 20 multiple array reads are 
averaged for each loop to reduce read noise and adjacent 

^ The CIRPASS data reduction software is available from 
|http : //www.ast . cam, ac .uk/^optics/cirpass/docs .html 



loop averages are used to identify and reject cosmic rays. 
Bias subtraction was also performed at this stage. Since the 
target was stepped across different lenslets in the IFU for 
each of the eight observations, a sky frame was created for 
each exposure from the seven other integrations and this 
was used to perform first order sky subtraction. The spectra 
were then extracted using an optimal extraction algorithm 
(Johnson, Dean & Parry 2002) to account for crosstalk be- 
tween adjacent spectra on the detector and to trace the cur- 
vature of the spectra from each fibre. The extracted spec- 
tra were fiat-fielded using a dome fiat, in order to remove 
the variation in throughput between fibres and the varia- 
tion in sensitivity between pixels. The dome flat had pre- 
viously been read- averaged, bias-subtracted and extracted. 
The data were wavelength calibrated using an argon lamp 
exposure, processed in the same way as the science data, 
again using optimal extraction. We used 17 argon lines and 
a quadratic fit to the dispersion, producing rms residuals of 
0.2 A (0.1 pix). Residual sky-lines were removed by fitting 
a low-order polynomial function to the wavelength-rectified 
data frames. Finally 3D (x,y,A) data cubes were constructed 
for each integration before combining the eight individual 
exposures using the known telescope offsets. We observed 
two standard stars Hip257 and Hipl06522 (J « 9 mag) to 
determine the flux calibration (these produced consistent 
calibrations). The observations of the standard stars were 
processed in the same way and used to flux calibrate the 
data cube. 



2.2 HST Photometry 

HST WFPC2 imaging was obtained from the archive (pro- 
gramme G05996, the CFRS survey). There were five I- 
band F814W exposures with a total exposure time of 6700 s. 
We re-reduced the data from the archive starting with the 
pipeline processed (fiat-fielded and dark-subtracted) indi- 
vidual exposures. These were then combined with integer 
pixel shifts determined from the world coordinate system 
and with cosmic ray rejection at the 3cr level. 

The GIM2D software package (Simard 1998; Simard et 
al. 1999) was used to determine an accurate inclination, 
fiux and scale length for the galactic disc in the HST 
F814W image. The approach of gim2d is to convolve ide- 
alized two-dimensional galaxy models with the instrumen- 
tal point spread function (PSF) and to subtract this from 
the original two-dimensional image: the best-fitting model is 
found by altering the model parameters to yield the smallest 
residuals from the real image. 

Synthetic instrumental PSFs were created by Ti nyTim 
v5.o (Krist 1995), which has the advantage that they are 
noiseless, can be sub-sampled and are created at the required 
location on the array. We produced PSFs sub-pixelated by a 
factor of 5 for WFPC 2 / WF 2 (F814W). We considered only 
an exponential disc profile (log / cx 1/r): as we are working 
in the rest-frame B-band at z = 0.8, the red bulge should 
produce negligible deviations from the fit of the exponential 
disc to the WFPC 2 F814W /-band data. We excluded the 
two H II regions which were bright in the rest-/? fsection l3.11 
from the gim2d fit to the surface brightness profile, by mask- 
ing regions of width 0.5 arcsec centred on the two knots. We 
experimented with different mask sizes and found this did 
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Figure 1. The upper panel shows the spectrum for knot A (soUd 
hne), the lower spectrum is for knot B (solid line), both for ex- 
traction apertures of ~ 1 arcsec X 1 arcsec. The centre of the Ha 
emission for knot A (solid vertical line) is clearly offset from the 
central wavelength for knot B (dashed line) . The dotted spectrum 
shows the sky emission. Sky subtraction residuals are still visible, 
although by selection the Ha line for this galaxy lies away from 
the sky lines. 



not significantly affect the disc parameters and these errors 
are included in the uncertainty in the rotation velocity. 

The galaxy centre, position angle on the sky, disc scale 
length, inclination and total flux were free parameters in the 
model. The best-fitting disc model from gim2d had a half- 
fight radius of 1.25 arcsec (10 kpc): the disc itself had a po- 
sition angle on the sky of 74 ± 2° , the inclination to the line 
of sight was measured to be 78 ±2° (sin i = 0.98). The ex- 
ponential disc scale length was 0.75 ±0.05 arcsec (5.7 ±0.4 
kpc). NICMOS NIC2 imaging for this galaxy also exists. 
We have examined this and find that the centre determined 
from the _ff-band imaging (which traces the older starlight) 
agrees with the centre presented here. 

Using SExtractor (Berlin & Arnouts 1996) the galaxy 
was found to have Ivoga = 21.7 with the H ii regions con- 
tributing ~ 25 per cent of the total flux. Galactic extinc- 
tion was determined using reddening E{B — V) =0.05 from 
Burstein & Heiles (1982). A correction of 0.88 mag was made 
for the inclination-dependent internal extinction, including 
the residual absorption of a face-on galaxy (0.27mag), using 
the prescription of TuUy & Fouque (1985) with an optical 



depth T =0.55. The magnitudes were adjusted by a small k- 
correction (a dimming of 0.02mag for a Sbc spiral template 
spectrum from Coleman, Wu & Weedman 1980) to match 
the observed /-band to the rest-frame B-band. 

3 DISCUSSION 

3.1 Distribution of star formation 

In this paper we present a map of the distribution of star- 
forming regions within this z = 0.819 galaxy. Unlike long- 
slit approaches this is a complete census of the spatially 
extended star formation. As can be seen from Fig. H the 
Ha emission is concentrated in just two star-forming regions 
and these are spatially coincident with the 'knots' seen in 
the HST /-band image (which corresponds to the rest frame 
B). From the HST image both of these 'knot'-like H ii re- 
gions are spatially resolved and noticeably elongated in the 
plane of the galactic disc. Knot A is near coincident with the 
projected centre of the galactic disc and extends over ~ 2.7 
X 2.4 kpc FWHM (after deconvolution with the WFPC2 
PSF). Knot B is about 2 arcsec from the centre and is ~ 
2.8 X 1.8 kpc FWHM (deconvolved). Their properties are 
tabulated in Table Q In determining the rest frame B-hand 
luminosity of the star-forming knots we have first subtracted 
off the best-fitting galactic disc model provided by gim2d 
(section 2.2). The surface brightness within the half-light 
radius of the star-forming knots in rest frame B correspond 
to 10^ Mq kpc~^ within the range given in Meurer et al. 
(1997) for starburst galaxies at < 2 < 3.5. We note that 
knot B (the H ii region with the largest projected separa- 
tion from the centre) lies exactly on the bolometric surface 
brightness versus angular frequency relation given in Meurer 
et al. (1997), based on the sample of Lehnert & Heckman 
(1996). We measure Ha at the greater than 5a level for each 
knot with fluxes of ~ 10~^^ ergs cm~'^ s~^ each (Fig.0. We 
note our total Ha luminosity of 2.4 x 10~^^ erg s~^ cm~^ is 
slightly less than the value of 4 x 10~^^ obtained by Glaze- 
brook et al. (2001) in their CGS4 UKIRT spectroscopy. In 
the absence of reddening our Ha fluxes would correspond 
to star formation rates of ~ 3Mq yr~^ per knot, where we 
have adopted the prescription of Kennicutt (1998), applica- 
ble for a Salpeter-like initial mass function (IMF). We note 
that the rest frame B fiux densities (Table would un- 
derestimate these star formation rates (for the same IMF) 
by a factor of ~ 2, which we attribute to the differential 
dust reddening between 4400 and 6563 A or to the fact that 
the B-hand flux, unlike Ha, is not a direct measurement of 
current star formation rate. We spectrally resolve the line 
widths within the individual H ii regions (Fig.[TJ. After de- 
convolution with the CIRPASS instrumental line width we 
measure o-„ of 44 ± 15 and 74 ± 10 km s~^. 



3.2 Large scale kinematics 

Using the Ha emission line, we measure a velocity separation 
Vsep = 220 ± 10 km s"^ for knots A and B (Table[Tll. How- 
ever, we do not measure a rotation curve for this galaxy: we 
simply have two discrete line emission regions. Coupling our 
CIRPASS IFU 3D spectrum with the HST images (specifi- 
cally the location of the H ii regions and the best fit to the 
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Figure 2. Left: the HST F814W image (corresponding to tiie rest frame _B-band at 2 = 0.8). Tlie field lias been trimmed to 6 X 11 
arcsec. The two prominent H II regions are clearly visible against the extended highly inclined galactic disc. Middle: the HST F814W 
image with the best-fitting disc model from GIM2D subtracted. Right: these regions are seen in our CIRPASS near-infrared spectrum. We 
show a 2D representation of our 3D data cube, collapsed along the spectral direction and showing the ~15 A around the Ha emission. 
The spatial scale and orientation are identical to that of the HST WFPC2 image. 



Table 1. Properties of star-forming regions 
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(9) 


knot A 


2.8±1.0 


24.11 


1.00 


9.4 


11935.9 ±0.4 


7.9 ±1.0 


40 ±20 


74 ±10 


2.7 


knot B 


12.6±1.0 


24.10 


1.01 


14.7 


11927.5 ±0.4 


5.8 ±1.0 


180 ±20 


44 ±15 


4.2 



The columns are as follows: (1) projected distance of star-forming knot from centre of galactic disc, (2) total magnitude (MAG_AUTO 
from SExtractor, Berlin &; Arnouts 1996) in model-subtracted 7-band image, (3) rest-B flux density from model-subtracted B-band 
image, (4) Ho line flux, (5) central wavelength of Ha emission, (6) FWHM of Ha line, (7) velocity shift with respect to centre, (8) 

velocity dispersion av (deconvolved), (9) star formation rate from Ha. 



centre of the galactic disc) we can set a robust lower limit of 
180 ±20 km s~^ to Vmax (where Vmax is defined as half the 
total velocity shift across the galaxy) . Of course as we see no 
evidence for a 'turnover' in the rotation curve Vmax could be 
significantly higher. We note that our lower limit on Vmax 
is significantly greater than the value of 120 ± 10 km s^^ 
presented by Barden et al. (2003) for the same galaxy based 
on their ISAAC/VLT near-IR spectroscopy. This is ~ \vsep 
and would be consistent with our measurement of the red- 
shift difference between the H 11 regions if it was assumed 
that the centre of rotation was the geometric midpoint of 
the two star-forming knots. However, the HST image does 
not support this, with knot B at far greater projected dis- 
tance from the disc centre. The systemic redshift is likely 
to be close to the Ha redshift of knot A [z = 0.8187 cf. 
z = 0.8191 from the CFRS catalogue, Lilly et al. 1995b). 

A simple estimate of the mass of the galaxy, M = 
V'^R/G, yields a value of 0.9 x lO" Mq within 12 kpc. 



3.3 The Tully-Fisher Relation at z ~ 1 

The TF relation is a fundamental scaling law, the evolution 
of which is critical to our understanding of mass assembly 
and star formation in disc galaxies. Evolution of the TF re- 
lation at 2 > 0.5 has been investigated with varied results. 
Vogt et al. (1996, 1997) find only moderate luminosity evo- 
lution (AM < 0.4) out to 2 ~ 1 while Barden et al. (2003), 
Milvang- Jensen et al. (2003) and Bohm et al. (2004) detect 
a brightening of ~ 1 mag in the rest-frame _B-band. Much 
stronger evolution is suggested by the study of Simard & 
Pritchet (1997), which finds disc brightening of ~ 2 mag at 
lower redshift {z ~ 0.35) and Rix et al. (1997), which derives 
a brightening of 1.5 mag a,t < z >= 0.25 using the [oii] line. 



-20 -21 -22 

M (mac) 



Figure 3. The rest frame B-band TF relation showing the po- 
sition of CFRS 22.1313 with B-band flux from the star-forming 
regions excluded (filled circle) and included (open circle), com- 
pared with high redshift (0.6 < 2 < 1) data from Bohm et al. 
2004 (crosses, high quality data only), Vogt et al. 1996 (stars), 
the field galaxy sample from Milvang-Jensen et al. 2003 (trian- 
gles) and the Barden et al. 2003 data point for CFRS 22.1313 
(diamond). IFU data for a 2 = 1 lensed arc from Swinbank et 
al. 2003 is also plotted (open square). The dashed line shows the 
high redshift TF relation of Barden et al. (2003). The dot-dash 
line shows the local TF relation of TuUy and Pierce (2000). 



The degree of evolution observed for a particular sample is 
likely to depend strongly on the selection method. 

In Fig. El we use our measured lower limit on Vmax to 
place CFRS 22.1313 on the TF relation and compare with 
the TF relation at the current epoch (TuUy & Pierce 2000). 
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We also plot the estimated TF relation at redshifts 0.6 < 
z <l{Az = 0.2 either side of the redshift of CFRS 22.1313) 
derived by other groups using long-slit spectroscopy. Using 
the total HST F814W magnitude for the galaxy of Jvcga ~ 
21.7 we obtain a lower limit on the offset of this one galaxy 
from the local rest B-band TF relation consistent with disc 
brightening of no more than ~1 mag a.t z = 0.8 and entirely 
consistent with no evolution. 



4 CONCLUSION 

Our main results may be summarised as follows: 

(i) We have presented the first near-IR integral field spec- 
trum of a 2; ~ 1 field galaxy. We have used our CIRPASS 
integral field spectrograph to map the spatial and velocity 
distribution of Ha in the z = 0.819 galaxy CFRS 22.1313. 

(ii) The Ha emission originates from two giant H ii re- 
gions offset by 220 ± 10 km s~^ : we do not observe a full 
rotation curve for this galaxy. 

(iii) Archival HST imaging shows a clear disc morphology, 
along with two bright star-forming knots coincident with the 
Ha emission and accounting for ~ 25 per cent of the total 
/-band flux (the rest frame _B-band). 

(iv) We have coupled the 3D spectroscopy with the best- 
fitting model to the underlying disc in the HST F814W im- 
age, and used this to derive a lower limit to the rotation 
velocity of 180 ±20 km s~^. The lower limit on the offset 
of this one galaxy from the local rest B-band TF relation is 
consistent with disc brightening of no more than ~1 mag at 
z = 0.8 and entirely consistent with no evolution. 

(v) Integral field spectroscopy eliminates many of the sys- 
tematic uncertainties inherent in long-slit studies. A larger 
sample of IFU data is required to accurately determine the 
evolution of the TF relation out to 2; ~ 1. 
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